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In this paper we describe how to derive the expressions for the higher nonlinear generation of waves, 
their transmission and reflection for the case of a normal incidence plane wave by direct superposition of 
dipoles oscillating anharmonically under the influence of the nonlinear polarization. We describe explicitly 

B 

that the transmitted nonlinear harmonic wave can only propagate inside the material if the phase matching 
condition is fulfilled and show for the case of second-harmonic-generation that our calculation yields similar 
results with coupled-mode-theory (CMT). Furthermore, we show using dipole superposition how to obtain 
nonlinear Fresnel equation for the case of reflection and discuss its implication if there is absorption . 
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I. INTRODUCTION 



Since the discovery of femtosecond lasers, nonlinear optics has become one of the standard 
tools for material analysis. Nonlinear optical laser techniques are attractive because they provide 
a high temporal, spatial, and spectral resolution, as well as its applicability to all interfaces ac- 
cessible by light 7 . The theory for nonlinear harmonics generation has been well understood both 
classically and quantum mechanically. 

The standard classical method to obtain expressions of transmission and reflection of the field 
inside a dielectric is usually performed by applying the macroscopic Maxwell equation on the 
boundary between two mediums where one will then obtain the Fresnel equation 7 7 7 . From the 
microscopic classical point of view, the dielectric consists of polarizable atoms or molecules, each 
of which is radiating in vacuum in response to the incident field and in response to the fields 
radiated by the other atoms. The total radiated field in linear optics is then the superposition of 
the incident field with all radiated fields. The incident field, however, drives at the same time, due 
to their polarisability, the microscopic dipoles. The superposition of their radiation fields, adding 
the incident field too, generates the correct reflected and transmitted field. Mathematically this is 
expressed as an integral equation for the field either inside or outside the material. The equivalence 
can be seen e.g. by comparing near field, far field, dispersion relation and Fresnel formulas for 
stratified layers, etc. 7 7 7 . 

II. EWALD-OSEEN EXTINCTION THEOREM 

Microscopically, nonlinear harmonics generation can be seen as radiation from the charges 
that experience anharmonic oscillation due to an external driving field. This view is actually an 
extension of the idea in linear optics that an external field impinging on a dielectric material pro- 
duces dipole oscillations in the material and their superposition constitute a radiating field. Using 
this approach, Ewald 7 and Oseen 7 were able to demonstrate linear refraction and transmission 
phenomena inside the dielectric microscopically, yielding an equivalent result if calculated from 
the solution of the macroscopic Maxwell equation. Although the macroscopic Maxwell equations 
yield the correct solutions, they hide the physics. 

For the linear case it has been shown by Fearn et. al. that a normal incidence plane wave 
electric field that propagates inside the material induces dipole oscillations that radiate waves at 
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Figure 1 : An incident plane wave from vacuum into a dielectric material is producing a field polarization 
along the x direction. This in turn produces dipole oscillation pararel to the polarization. 

the speed of light, c inside the material. The integral equation describing the total field at a depth 
z from the surface derived by Fearn et. al. ? reads. 



oo 








In linear optics this is an integral equation, which has to be solved. For the case of nonlinear 
harmonic generation, however, the mathematical structure of the underlying equation to be derived 
will be, as seen below, much simpler. This is due to the following facts: a) There is no incoming 
nonlinear incident field, all the nonlinear fields are generated by the anharmonic oscillating dipoles 
driven by the nonlinear polarization inside the material and b) the amplitude of the nonlinear field is 
far smaller than the amplitude of the driving one. Strictly speaking, if the nonlinear field amplitude 
obtains the same size than the driving one, coupled integral equations between the linear and the 
nonlinear field(s) would result. However, in the low depletion approximation we assume that the 
amplitude of the fundamental field is space independent, and then, instead of an integral equation, 
only an integration has to be performed. Surprisingly enough, to the best of our knowledge no one 
has calculated this nonlinear dipole superposition explicitly even for the normal incidence case. 

We proceed now in the following way. We first derive, taking the approach of Fearn 7 , the 
equation for the nonlinear field. Then we derive for a material with finite thickness L, the su- 
perposition of the nonlinear radiation fields at + and L and compare it with results obtained 
from coupled-mode-theory (CMT). The reflection due to integration of all nonlinear dipoles is 
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also calculated. 



III. SUMMING OVER NONLINEAR DIPOLES 

Here we present a way to obtain an expression for the transmitted and reflected nonlinear field 
that is produced by superposition of all radiating (nonlinear) dipoles inside the dielectric. In doing 
so, we follow a similar approach as in Ref . 5 but use the nonlinear polarization instead of the linear 
case. It has to be noted however that one should not expect to obtain the same results as in the 
linear case because no second harmonic field is incidence on the material and that the nonlinear 
dipoles are driven by the nonlinear polarization. For simplicity, let us assume the external field is 
sufficient large to generate SHG inside a dielectric and that it is incoming at normal incidence. 

To obtain the nonlinear field we calculate the contribution from all nonlinear dipoles. We 
start by expressing the electric field produced by a single dipole (Hertz dipole) at r' inside the 
dielectric 7 : 

R iP \ C m J ^ C m R 2 dt 



E d (r,t) = [3{p.n)n-p] 



+ 

(2) 



r R dt 2 

Here p is the unit vector for the dipole polarization, R = \r — r'\, and h is the unit vector directed 
towards r — r'. Note that in contrary with the linear case, the retarded term which 
describes the time evolution of the dipole field is propagating at the speed c rn = c/ y/e(u) instead 
of c, because the dipoles are driven by the nonlinear polarization that is also propagating with 
c/ y/e(u) inside the material. 

The linear case point of view can be extended to cope with nonlinearity by considering that if 
the field intensity is sufficiently high, the nonlinear far field emerges from the anharmonic motion 
of the charges affected by the nonlinear driving force of the incomming electric field. As in 
the linear case a nonlinear polarization results from this process. The nonlinear polarization can 
be obtained clasically by expanding the harmonic oscillator Lorentz model assuming nonlinear 
restoring forces excerted on the electron 7 : 

r + 2^r + ujir + (3 2 r + (3 3 r + ■■■ = — (3) 

m 

where to, e, r are respectively the electron mass, charge, and its displacement from equilibrium, 
— I^/ttlt is the damping force and the restoring force —uolmr — l3 2 mr — f3 3 m'r + ... describes 
the nonlinear effects (second harmonic, third harmonic, ect). Because the linear field E w is far 



higher than the nonlinear fields we can neglect the higher nonlinear driving fields. Here (3 is the 
hyperpolarizability (tensor). The complete expression for the polarization consist therefore of the 
linear and nonlinear parts: 

p = ^E[ oc {z')e- iwt + p 2 E[ oc {z'fe' i2wt + hE\ oc {z'fe^ ut + ... (4) 

or in a more compact form: 

p = J2PnE l r(zTe- in0Jt (5) 

n=l 

where n denotes the harmonic order. Note that (3 is aside of the linear part (first term of eq. (4)) a 
tensor and that we write the nonlinear local field here as a multiplication factor, m of the external 
driving fundamental field so we set E[ oc (z') n = mE 1 (z') n . We next calculate the total (linear and 
nonlinear) dipole radiation due an incoming field at normal incidence (FIG. 1). Inserting eq (4) 
in (2) and by noting that p = x and h ■ x = (x — x')/R then integrating over all volumes and 
summing over all the terms yields: 

E d (r,t) = J2mNP n fffdx'dy'dz' [E^z'T e~ ik ™ R \(?> { ^ - l) - - 

LV J (6) 
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Here we have used the notation for the higher harmonic wave vector k nuJ = nk. It turns out 
that the integrals in eq. (6) is easier to solve in polar coordinates by introducing the coordinate 
transformation x — x' = pcos(p and p 2 = R 2 — (z — z') 2 . The mathematical methods to perform 
the integration follows directly from 7 but by using a generalized wave vector k nu] instead of k : 

POD 

E d (z) = Y^2mmk nuJ N(3 n / dz 1 x [E 1 (z')] n e ik ^ z ~ z '\ (7) 

n=l J ° 

Eq. (7) is the total dipole contribution from the linear and nonlinear harmonics and can be used 
to calculate the transmission and reflection formulas for nonlinear optics. Note that there is a 
harmonic propagation given by e %k ^\ z ~ z '\ which is similar to the results using antenna dipole 
approach in calculating second harmonic generation (SHG) ? . For SHG we simply take the second 
term in the summation and set n = 2 in eq. 7: 

PCX) 

E d2 M) = 2mmk 2uJ N(3 2 / dz 1 x [E^z')} 2 e lk ^ z ~ z '\ (8) 

Jo 

Before we proceed further we stop to state an important assumption which will significantly sim- 
plify the calculation. For many cases, such as studying SHG in centrosymmetric media, the non- 
linear field E 2 is far smaller than the linear field Ex. Therefore, the calculation can be performed 



independently (e.g. it does not require self consistency) for each nonlinear polarization term. This 
assumtion is similar to the low depletion approximationt which allows couple mode calculations 
be analyzed independently (uncoupled). 7 7 

IV. NONLINEAR TRANSMISSION 

To calculate the nonlinear transmission for the n-th harmonics we integrate eq. (7) over z > 
and substitute [Ei(z')] n = [£' 1 e lfez '] n = E™e mkz ' where k nuJ is the n-th harmonic propagation 
vector: 



The integration above can be better understood by looking at Fig 2. The field at a point z' is 
obtained by superposition of all dipoles from the left which is represented by the first integral in 
Eq. 9 and superposition from all dipoles from the right represented by the second integral. It is 
also worth to notice that the second term in Eq. 9 is very small because k and k nuJ is in the order 
of 10 7 indicating that the superposition from all dipoles from the right can be neglected hence 
the contributions will come from the dipoles to the left. Also, in order for a nonlinear wave to 
propagate inside the material requires AA; = k nu] — nk must be small or in other words the phase 
matching condition is necessary. 

We now consider the case where a nonlinear incomming normal plane wave passes through a 
material with a coherence length L (FIG. 2). Calculation of the nonlinear field is usually solved 
either by coupled mode theory 7 or by means of the dipole antenna method 7 . Because from Eq. 9 
we know that the contributing dipoles comes from the left we need only to integrate from to L. 
Hence the field after propagating through the coherence length L produces a transmission: 



The part in the pharanthesis denotes the evolution of the wave envelope whereas the remaining 
describe the phase propagation of the harmonic wave. The intensity is simply obtained by squaring 
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(9) 




(-2nmk nuJ N(3 n E™ [e k ™ z/2 Sinc(AkL/2)]) e 



(10) 



dz' 




Figure 2: The nonlinear far field transmitance for a normal incidence fundamental field can be calculated 
by summing all the radiation from the anharmonic dipoles affected over a coherence length L. 



the field envelope 

I nuJ - oo^Sinc 2 (AJfeL/2) 
for the case of second-harmonic-generation (SHG) we set n = 2 in eq. (10) 



/*„ - u 2 2 ElSinc 2 {AkL/2) 



<-2w 
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(12) 



which is the celebrated intensity formula obtained from coupled mode theory (CMT) or the an- 
tenna picture. 

V. NONLINEAR REFLECTION 



For the reflection eq. (7) has to be analyzed with z < or taking the second integration in eq. 
(9) and evaluate from to infinity 



E d (z) nuJ = 2irimk w NP n / dz' x E ^e ikz ' e ik ^ z '^ 

Jo 

-27tmk nuJ N /3 n E? 



(13) 



j\k -\- k nul f 

Interestingly, we can derive the nonlinear Fresnel formula by aplying the Lorentz-Lorenz relation. 
Using this the local field for second harmonics (n = 2) can be expressed as an average of the 
driving field by a correction of: 



Nf3 2 E l ° c = -L [e(2u) - 1] E x 

47T 



(14) 



therefore for SHG we put the multiplication factor m = 47r ^- fe \n(2oo) 2 — 1] and set n = 2 in eq 
(15) 



E d {z) 2uJ = ^ [n{2ujf - I] 2 J dz' x E X E\ 



loc^i2kz' gik2tu (z'—z) 



[n{2u) 2 - 1] Joc f /■••>., \ ,■/,,. : 



2 ^UtU^ (15) 

in reflection, which is measured in vacuum (air) there is no dispersion, thus we can set k 2u = 2k 
or in other words there is always phase matching, thus we can write 

2k 1 



2k + 2k [n (2u) + 1] 
and using the relation given in Ref. 5 



(16) 



Ei = 2 ^^E i = — -Ei (17) 

1 l + n(w) l + n(2w) 

where is the incidence fundamental field. Inserting eq. (16) and eq. (17) in eq. (15) we obtain 
the Fresnel Formula for SHG in reflection 

E d (z) 2uJ = n ^ U ] - 1 E t E\ oc e- ik2 " z (18) 
dK J2uJ n(2w) + l 1 

Similar result can also be obtained for higher orders (n = 3, 4, etc.). We then proceed by dis- 
cussing the presence of absorption. It is well known that in absorbing medium the driven and 
nonlinear wave experiences an exponential decay which can be incorporated in our contex as 

[Ei(z')] n = [Eie tkz ' e~ nKZ '] n = [E 1 e mkz ' e~ KZ '] and e tk «^( z '~ z ) e ~ ,i2Z ' respectively, where k is the 
extinction coefficient. In this case the reflected nonlinear field inside the dielectric becomes: 



E d (z) nuJ = -2nmk nui Np n E? 



(19) 



where k s = nn + k 2 is the total extinction due to absorption. Note that if absorption is small 
< k s < 1 the reflection will increase whereas if the absorption is very high (k s >> 1) the 
reflection will decrease significantly. 



VI. SUMMARY 



We show by direct superposition of dipoles that the formulas of reflection and transmission 
for the nonlinear harmonics can be obtained. The results also show that phase matching is an 
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important requirement for nonlinear wave propagation and that we can obtain the nonlinear Fresnel 
equation for the case of reflection as well as deduce implications of absorption. 
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